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Abstract 

In order to perform numerical simulations of the KPZ equation, in any dimensionality, a spatial 
discretization scheme must be prescribed. The known fact that the KPZ equation can be obtained 
as a result of a Hopf-Cole transformation applied to a diffusion equation (with multiplicative noise) 
is shown here to strongly restrict the arbitrariness in the choice of spatial discretization schemes. On 
one hand, the discretization prescriptions for the Laplacian and the nonlinear (KPZ) term cannot 
be independently chosen. On the other hand, since the discretization is an operation performed 
on space and the Hopf-Cole transformation is local both in space and time, the former should be 
the same regardless of the field to which it is applied. It is shown that whereas some discretization 
schemes pass both consistency tests, known examples in the literature do not. The requirement 
of consistency for the discretization of Lyapunov functionals is argued to be a natural and safe 
starting point in choosing spatial discretization schemes. We also analyze the relation between 
real-space and pseudo-spectral discrete representations. In addition we discuss the relevance of the 
Galilean invariance violation in these consistent discretization schemes, and the alleged conflict of 
standard discretization with the fluctuation-dissipation theorem, peculiar of ID. 
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I. INTRODUCTION 



Soon after its formulation in 1986, the KPZ equation 



A 



m 



dth = iydih + - {d^hy + F + e^{x, t), 



(1) 



became a paradigm as the description of a vast class of nonequilibrium phenomena by means 
of stochastic fields. The field h{x, t) whose evolution is governed by this stochastic nonlinear 
partial differential equation, describes the height of a fluctuating interface in the context of 
surface growth processes in which it was originally formulated. In particular, Eq. ([T]) assumes 
a one-dimensional (ID) homogeneous substrate of size L. The parameter u determines the 
surface tension, A is proportional to the average growth velocity (the surface slope is parallel- 
transported in the growth process), and F is an external driving force. Finally, ^{x,t) is a 
Gaussian white noise with {^{x,t)) = and {^{x,t)^{x' ,t')) = 26{x — x')6{t — t'). As usual, 
periodic boundary conditions are assumed. 

From a theoretical point of view the KPZ equation has many interesting properties, 
like its close relationship with the Burgers equation j4| or with a diffusion equation with 
multiplicative noise, whose field 0(x, t) can be interpreted as the restricted partition function 
of the directed polymer problem. But clearly, investigating the behavior of its solutions to 
obtain e.g. the critical exponents in one or more spatial dimensions js-jlS] requires the 
(stochastic) numerical integration of a discrete version. Although a pseudo-spectral spatial 



discretization scheme has been recently put forward 14l4l6l|. as well as a numerical large 



deviation theor y ^ \ , real-space discrete versions of Eq. ([T]) are still largely used for numerical 
201 ]. because of their relative ease of implementation and of interpretation 



simulations 



17 



in the case of non-homogeneous substrates (for instance, a quenched impurity distribution 
2l| ) among other reasons. To that end, several real-space discretization schemes have been 



proposed [10 
simulation. 



y 



19| . which are claimed to cure particular "diseases" of the numerical 



In the present work, no attempt is made of comparing alternative real-space discretiza- 
tion schemes in sought of special KPZ features. Instead, we seek to point out some basic 
conditions that any spatial discretization must fulfill in order to consistently describe the 
KPZ equation. Nonetheless, for the sake of brevity and for ease of comparison with other 
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proposals, we shall adopt the notation in Ref. 19|, namely (calling Ax = a) 



L = aN, ^ hHi^h^_ 



and 



2(1 + 7) 

with 7 G [0, 1]. On one hand, the restriction to A;, / G {—1, 0, 1} is unnecessary. On the other 
hand, we shall denote Lj — )■ Ln\ (hj) and A^j^^ N^^\hj). [The subscript (1) indicates 
that only nearest neighbors are involved in the prescription of the discrete Laplacian.] By 
analogy, we shall write 

*Sfc = i^j+i - 4>j+k)/a, 



and consequently 



and 



iVW(0,) = 



2(1 + 7) 

Two main symmetries are usually ascribed to the ID KPZ equation: Galilean invariance 
and the fluctuation-dissipation relation. 

• The flrst one has been traditionally linked to the exactness (in any spatial dimension- 
ality) of the relation a + z = 2 among the roughness a and dynamic z exponents 

y 



23|, although this interpretation has been recently criticized in other nonequilib- 



rium models 



25| . The roughness exponent a characterizes the surface morphology 
in the stationary regime t ^ t^. On the other hand, the correlation length scales 
as ^(t) ~ t^/^ with the dynamic exponent z, and is the time at which it satu- 
rates, namely C,{t ^ t^) ~ L. The ratio /3 = a/z is called "growth exponent" and 
characterizes the short-time behavior of the interface. 

The second symmetry essentially tells us that in ID, the nonlinear (KPZ) term is not 
operative at long times or in other words, that the long-time ID interface is equivalent 
to a path of Brownian motion jsl. From a theorem by Kolgomorov, this implies that 
the interface is Holder continuous with exponent strictly smaller than 1/2. For higher 



dimensions, the KPZ roughness exponent a decreases, implying a loss of regularity. 
Hence the error terms of a local numerical method (as e.g. a finite differences scheme), 
which are proportional to some higher-order derivative of the field, are not controlled. 
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As a consequence, a global method such as a pseudo-spectral scheme 
more adequate. Nevertheless, previous experiences found in the literature showed 
that finite differences schemes are still able to capture the universal features of KPZ 
evolution. This, together with our previous considerations, is our motivation for the 
present work. 

In Sec. [TT]we show that the relationship established by the Hopf-Cole transformation — 
between the KPZ equation and a diffusion equation with multiplicative noise [2I] — poses 
constraints on the discretization procedure. We verify the consistency of the standard 
(nearest-neighbor) discretization scheme and find the form of the corresponding KPZ term 
for a general real-space discrete Laplacian; we also present some comments regarding the 
mapping of KPZ into the directed polymer problem 2|]. In Sec. IIIII we analyze the problem 
from the perspective of the Lyapunov functional, show in what sense known prescriptions 
for the KPZ term from the literature fail the test, and find the corresponding consistent 
prescriptions. Moreover, we propose a consistent real-space discretization scheme whose 
accuracy is far higher than that of schemes of similar complexity in the literature. In Sec. 
IIVI we discuss the relation with the pseudo-spectral method. In Sec. |Vl we show that a 
consistent discretization scheme does not (essentially) violate the fluctuation-dissipation 
relation, peculiar of ID, and discuss the role of the Galilean invariance for the discrete rep- 
resentations of the KPZ equation, showing that such invariance seems not to be a necessary 
element to define the KPZ universality class. In Sec. I VI I we discuss a recently introduced 
variational approach for the KPZ equation {26 1, and show that it offers a natural framework 
for its consistent discretization. In Section IVIII we present some numerical results regarding 
critical exponents and the violation of Galilean invariance. Section IVIIII contains the con- 
clusions and final discussions. It is worth here commenting that some preliminary results 



were presented in 



27| 
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II. THE LAPLACIAN DETERMINES THE NONLINEAR TERM 



In this section we elucidate — by considering the standard, nearest-neighbor discretization 
prescription as a benchmark — one of two constraints to be obeyed by any spatial discretiza- 
tion scheme. It is very important to remark that this constraint arises due to the mapping 
between the KPZ and the diffusion equation (with multiplicative noise) through the Hopf- 
Cole transformation. Hence, for a general real-space discrete Laplacian, we state the form 
of its corresponding KPZ term. Even though the present analysis is performed on the KPZ 
equation, it is general in the sense that for sets of equations related among themselves 
through a local transformation there should be a consistent relation between the discrete 
transformed forms. 



A. The simplest case 



As it is known, the diffusion equation with multiplicative noise 



is related to the KPZ equation [Eq. ([T])] through the Hopf-Cole transformation 

X 



{x,t) = exp 



2z/ 



h{x, t) 



(2) 



(3) 



Note that this transformation isjust one particular example of the general implicit trans- 
formation written down in Ref. 



The standard spatial discrete version of Eq. ([2]), after transforming to a co-moving refer- 
ence frame (p —> (p + Ft, is 

Xe 

4- = z/L(i)(0j) + — 0^-, (4) 

with 1 < j < = 0, because periodic boundary conditions are assumed as usual (the 
implicit sum convention is not meant in any of the discrete expressions). The discrete noise 
C,j{t) is a Gaussian random variable with zero mean and correlation given by 



(e,(t)e.(t)) = 2^5(t-tO- 



Then, using the discrete version of Eq. (jS]) 

0j(t) = exp 



2u 



(5) 



(6) 
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we get 



namely 



Q2u 



Aa2 



31 



with S'^ = ^ifj^^. It is worth commenting here that this last expression was also pointed 
out in [9|, discussing aspects of discretization instabilities and the relation to the directed 
polymer problem. We will further discuss the mapping to the directed polymer problem in 
Sec. Ill CI below. By expanding the exponentials up to terms of order of a^, and collecting 
equal powers of a (observe that the zero-order contribution vanishes) we retrieve 



with 



Q(i){hj 



1 



■j+i 
3 



H 



(7) 



(8) 



{Q stands for "quadratic"). As we see, the first and second terms on the r.h.s. of Eq. ([7]) 
are necessarily related by virtue of Eq. ([6]). 



B. The general case 



A Taylor expansion of around shows that the general form of the discrete Lapla- 
cian, involving up to the n-th nearest neighbors of site j, is of the form 



3-1 



where as before, the subscript stands for the number of nearest neighbors. Since the maxi- 
mum value for n is M = (A^ — 1)\2, where \ denotes integer division, one may alternatively 
run the sum up to M and set 6; = 0, Z = + 1 . . .M. The remaining hi, that are other- 
wise arbitrary, should be fixed by whatever criterion (below, we shall use the criterion of 
maximizing accuracy). 

Repeating the steps described above, one obtains 



L{n){hj) 
Q{n){hj) 



3+1 



H 



H: 



-3-1 



(10) 
(11) 
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C. Few remarks on the directed polymer problem 



We devote this subsection to briefly comment about the mapping of KPZ onto the directed 
polymer problem. Such a mapping can be carried out via the Hopf-Cole transformation ^ 
and the resulting linear equation corresponds to Eq. In order to employ the usual rules 
of calculus, here we assume the Stratonovich interpretation for the multiplicative noise. The 
corresponding flnite differences scheme is, explicitly, 

XF Xe 

4>3 = ^^{i)(0i) + + (12) 

As indicated in Eq. ([5]), the discrete noise E,j{t) is a Gaussian random variable. The mean 
value of Eq. ([12]) is 

dt 2v ■' Aua 

One immediately realizes that the drift of this equation becomes singular in the continuum 



limit a — 7- 0, so one has to renormalize this theory |29H31|. This is done by decomposing the 
bare parameter into an effective and a singular component, F = F^s+Fs, with Fg = — l/(2a). 
The resulting equation is then 

dt ~ ^ 2v ^^^^ 

which is flnite, but in which F^^ has to be measured directly from the experiment. Thus the 
correct interpretation of the Stratonovich Eq. (IT^ is the following Ito equation 

bj+i + - 20J AFeff Xe 



V- 



+ — + (15) 



a2 2v Iv 



In order to measure the effective growth rate, one can solve the linear Eq. ( iMl) to flnd 
the globally stable solution {4>j(t)) = (0j(O)) exp[AFofj t/(2z/)] for a spatially homogeneous 
initial condition (0^(0)) = (0(0)). And so this effective rate can be measured from experi- 
mental/numerical data in the following fashion 

Feff = ^ln[(0,(t))/(0,(O))], (16) 

or alternatively 

i^efr = ^ln{(exp[A/i,(t)/(2i.)])}, (17) 
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assuming that the initial condition is hj{0) = 0. Applying Jensen's inequality to this last 
relation one finds 

{hjit))>F,st, (18) 
in agreement with what one could directly obtain from the KPZ equation [Eq. ([T])]. 



III. EXPLOITING THE DETERMINISTIC LYAPUNOV FUNCTIONAL 

An important feature of the Hopf-Cole transformation — Eq. or (jH]) — is that it is local, 
i.e. it involves neither spatial nor temporal transformations. Some effects of this feature are 
the following 

1. The discrete form of the Laplacian — namely the operator L(„) — is the same, regardless 
of whether it is applied to or to h. 

2. For a given (i.e. a given set of 6^ = 0, / = l...n), Q(^n) is determined by the 
Hopf-Cole transformation, Eq. ([6]). 

In this section, we want to go further with the criterion that the definitions of the discrete 
operators should not depend on the fields on which they are applied. 

The deterministic part of Eq. namely the diffusion term, admits a local Lyapunov 
functional. In other words, for e = 0, Eq. ([2]) can be written in the following variational 
form 

with 

^[0] = ^|da; (9.0)2. (20) 

The aforementioned criterion dictates the following set of discrete forms (thus Lyapunov 
functions, for any finite N) of Eq. (j20|) 

1 ^ 

J'inm = :,^aJ2Qin){<l>,). (21) 



It is a trivial task to verify that 



1 gJ^(n)[0] , . 

^I^(n)(0,)---^^. (22) 
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There is no loss of generality in taking j = = 0. If we rearrange the sum in Eq. (12T|) 
as T,jLM+i-N^ with M = {N - 1)\2, then only -n < j < n will contribute in Eq. fl22|l . 
Moreover, their contribution is such that they cancel the factor 1/2 in front of the sum in 
Eq. T2T\i . For completeness, let us show the particular functional form for Q{i){4>j) 



7/ ^ 



+ 



(23) 



Other discrete forms of the KPZ term 



Of course, Eq. fl22|) does not uniquely determine the Lyapunov function. Expressions 
other than Eq. (!2T|) may yield L(^n){4'j), provided that they contain the right terms, in the 
right proportion. Take as an example the proposal of Refs 
19| 



10 



Hi, coded as A^f^^^ in Ref. 



3 



i+1 



1 



(24) 



By using N^^/'^\(j)ji) instead of Q{n){4>j') in Eq. ( 120|) . we obtain 



1 



,i-2 



2^(i)(0.) + -i'^i 



i+2 



(25) 



which is an instance of L(2)(0j), with bi = 2, 62 = 1- As it was shown before, the procedure 

(1/ 

)(l/2)('/,A— Jo rH-J + J^A^ I (U:i-^\^ I I I ^26) 



outlined in Sec. HTl will yield L^^^'^\hj), together with 



and noi N^^^'^^hj). Hence, the proposal of Refs. 10|, lul is not consistent: On one hand, 
N^^/'^\hj) does not correspond with L(i)(/ij), as it is used. On the other hand, it does not 
correspond with L^^^'^\hj) either, as showm 



As stated before, the proposal of Refs. [lO|, lUl] belongs to a family coded in Ref. [l9|] as 



(7) 



7 e [0,1]. The choice 7 = 1 yields A^^^^^i) = i ("^^j-i)^- the other hand, the 
choice 7 = yields L^^^{hj) = L{i){hj). They all correspond to n = 2, with 61 = 1 — 7, 
62 = 7/2. The equivalent of Eq. (l24l) is now 

1 



AtW (</,,.) 



2(l + 7)a2 



(27) 
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which V7 G [0, 1] yields 



and 



1 + 7 
1 

2(1 + 7) 



:i - 7) {//f ^ + Hf') + I + H^-' 



{2L(i)(/ij) + 7 + } 



1-7) m^'r+{Hf') 



+ 
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However, the accuracy of this discretization is unknown and should be studied. 



B. A more accurate discretization scheme 



(28) 



(29) 



Again, a Taylor expansion of (pj+i around (pj shows that the C(a^) corrections to 
[applied to hj in Eq. ffTOl) ] are of the form 

Thus, the C(a^) correction to L^g] is j2^i+^^^- attains its minimum value (^c?f/i) 
precisely for 7 = 0, namely for L(i). What is then the convenience of a more complex 
prescription for the Laplacian? 

A wise criterion for choosing bi and 62 in -^{2) is making the O(a^) corrections vanish. 
This yields the prescription bi = 16, 62 = —1, known to be accurate up to corrections of 
Oia^) [321. 

Carrying out the procedure sketched in Sec. Ill Al we obtain 



Q{2){hj) 



12 



24 



(30) 
(31) 



The 0{a?) corrections to Q^n) are 



?>{dlhf + A{dlh) {d,h) 



4! Er=i I'bi 
which also vanishes for bi = 16, 62 = ~1- 

Since this discretization scheme fulfills the consistency conditions, is accurate up to 0{a^) 
corrections, and its prescription is not more complex than the ones studied before, it is 
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obvious that it will be a convenient one to be used when a higher accuracy in numerical 
schemes . .e,ui.ed. The pcibUity that tt may also help to conttol |ot at least delay) 

the numerical instabilities found in previous works (see for instance |9|, |33| and references 
therein) will be the subject of further work. 

Let us remind again that, as we already pointed out in the introduction, these results are 
formal since the higher order derivatives of the field are not under control. 



IV. PSEUDO-SPECTRAL DISCRETIZATION 



As was indicated in the introduction, a pseudo-spectral spatial discretization scheme has 



been recently introduced 



14-16 



In this section we show the relation existing between 
the present analysis and the indicated pseudo-spectral spatial discretization scheme. 

The pseudo-spectral discretization procedure starts by Fourier expanding the field h{x, t) 



°° - / 27r 

h{x,t) = hk{t)exp \A—kx] , 

k=—oo 



with 



1 



hk{t) = — dx h{x,t) exp —i—kx 
L Jo \ L 



,27r. 



and Eq. ([3]) can be rewritten as 

0(x, t) = exp 



^ Y hkit) exp (i'^kx 

k=—oo 



(32) 



(33) 



(34) 



Thus Eq. (ED reads 



Y exp [i^kxj I hkit) + i^^j khk{t) 



A;=— oo 



27r\- 



A °° / 27r 

vk -\ — ^ k'hk'{t)exp {'i—k'x 



k'=—oo 



-£a(t)j = 0, 

(35) 

since C,{x,t) is also assumed to be L-periodic as a function of x. A sufficient condition is 
that 



2tt 



hk{t) = - ( Y ) ^^'^(^^ 



vk + — Y, k'hk'{t) exp fi — k'x 



+eUt). (36) 

k'=—oo 

In this context "discretize" means to consider only Fourier modes, including k = 0. If 
M ={N - 1)\2, then 



hk{t) 



khk{t) 



2'jx 

vk + - Y k'hk'it) exp ( i—k' 

^ k'=M+l-N 



X 



+ eUt)- (37) 
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As indicated before, it is interesting to connect real-space and pseudo-spectral discretiza- 
tion approaches. From Eq. f l32|) [with h{ja,t) = hj{t)] we have 



H: 



j+i 



2 ^ ^ 

hk(t) sm 

^ k=M+l-N 



Txkl 



exp 



and 



M 



L{M){hj) = ^k{t) 
k=M+l-N 

By equating this expression to 

M 

- E Mt) 



cos 



f 2'Kkl \ 
\ N ) 



.27r, . 
exp ( i—kj 



(38) 



(39) 



27r . 
exp \ i—kj 



k=M+l-N \ / 

we might think of the pseudo-spectral discretization as a particular real-space discretization, 
whose coefficients are the solutions of the linear system 



M 
1=1 



2'Kkl 



1 (2T^kV 



This equation is linear and homogeneous, and so it admits the trivial solution bi = 0\f I & 
{1, ■ ■ ■ , M}. This equation expresses in fact the fundamental difference of the spectral and 
finite differences discretization: the lattice spectrum. For the finite difference scheme the 
discrete Laplacian is no longer (27rA;/L)^ (namely twice the spectrum in the continuum) but 
For X = L/2, already for k = 1 the difference is tt'^/2 — 2 fti 3. If we equate 



1 



cos 



2nkx 



instead Eq. ( 139|) to 



M 



E h,it) 

k=M+l^N 



COS 



27ckV 



- 1 



exp ( i—kj 



then Eq. (j^Oj) says nothing new. There is still a complete arbitrariness in the choice of the 
coefficients 6/. This corresponds to the fact that this scheme is nothing but the Fourier 
transformed version of the finite differences one. 

V. GALILEAN INVARIANCE AND FLUCTUATION DISSIPATION RELATION 



There are two main symmetries associated with the ID KPZ equation: the fiuctuation- 
dissipation relation and Galilean invariance. On one hand the fiuctuation-dissipation rela- 
tion essentially tells us that the nonlinearity is asymptotically (that is, at long times) not 



operative in ID. On the other hand, Galilean invariance has been traditionally re 
the exact relation among exponents a + z = 2, that holds for all spatial dimensions 
However, it is worth remarking that this interpretation has been recently criticized 



ated to 



22 



24|. 



23|. 



12 



A. Galilean Invariance 



Galilean invariance means that the KPZ equation is invariant under the transformation 



X — )■ X — Xvt, 
h —7- h + V X, 



(41) 



where v is an arbitrary constant vector field and F is the external (constant) driving force. 
Using the classical discretization 

d^h^^Hj^l (42) 
for the (complete) KPZ equation, we find 



A 



(43) 



One can immediately check that this equation is invariant under the discrete Galilean trans- 
formation 



j a j a — Xvt 
hj — )■ hj + V j a 
F ^ F 



X 



(44) 



However, Eq. ( H3|) has been criticized for its instability properties, at least when the spatial 
discretization is not fine enough {q]. If we use the alternative discretization 

A 



+ F + i^{t), 



(45) 



we find that this equation is not invariant under the discrete Galilean transformation. In 
fact, the transformation h ^ h + vja yields an excess term which is compatible with the 
gradient discretization in Eq. (142|) : however, this discretization does not allow to recover 
the quadratic term in Eq. fj45|) . indicating that this finite differences scheme does not fulfill 
Galilean invariance. The Hopf-Cole transformed equation 



XF X 



(46) 
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is Galilean invariant, i.e., it is invariant under the transformation indicated in Eqs. 
Hence, the nonlinear Hopf-Cole transformation is responsible for the loss of Galilean invari- 
ance. Note that these results are independent of whether we consider this discretization 
scheme or a more accurate one. 

Galilean invariance has been always associated with the exactness of the ID KPZ ex- 
ponents, and with a relation that connects the critical exponents in higher dimensions. If 
the numerical solution obtained from a finite differences scheme as Eq. fH5|) . which is not 
Galilean invariant, yields the well known critical exponents, that would strongly suggest that 
Galilean invariance is not a fundamental symmetry as usually considered. 

In Sec. IVIII we present some numerical results for the critical exponents using the con- 
sistent discretization schemes indicated in Eqs. ([H]) and flHTl) . and compare with those found 
with the standard one. All the cases exhibit the same critical exponents. Moreover, let us 
note that the discretization used in Refs. 



10 



ll| . which also violates Galilean invariance. 



yields the same critical exponents too. 

When we compare the classical discretization given by Eq. (H3l) . that explicitly reads 

^^-"^ a- +2 2a + 



(47) 



with the alternative one in Eq. ( H5|) . that reads 
h 



hi 



X 



H+i + hj^i — 2hj 

^5 + - 

a2 4 



h 



■j+i 



hj 



+ 



+ F + i^{t), 



(4^ 



we find that this second one presents excess fluctuations with respect to the first. This can 
be easily seen by means of the inequality 



{hj+i - hj.iY = {hj+i - hj + hj - hj^if < 2 {hj+i - hjf + 2 {h^ - hj^if 



(49) 



which immediately translates into 

A / hj+i - hj^i V 
2 I 2a 



X 

< - 
- 4 



hj+i - hj 



+ 



(50) 



where the inequality is strict unless hj = (/ij+i + hj^i)/2, an event which happens with zero 
probability (note that in ID and for long times, the KPZ interface has independent Gaussian 
distributed increments, as Brownian motion). This implies that the excess fluctuations are 
genuinely present in the interface dynamics. 
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The excess fluctuations from Eq. (HHj) respect to Eq. (H7|) can be explicitly computed: the 
alternative discretization scheme may be written as 



hj+i + h 



= ^— ^ + - \ — — - — - — !■ + 



A 



2a 



(51) 



where the term between curly brackets denotes the Galilean invariant fluctuations and the 
term between square brackets denotes the excess fluctuations. If the excess fluctuations 
are comparable to the Galilean fluctuations then there will be a strong violation of Galilean 
invariance. If the critical exponents still persist in this case, that would indicate that Galilean 
invariance is not such a fundamental symmetry as usually considered. This will be discussed 
in Sec. IVTH 



B. Fluctuation dissipation relation: stationary probability distribution 

As we have already mentioned, together with Galilean invariance, the fluctuation- 
dissipation relation is another fundamental symmetry of the ID KPZ equation. It is clear 
that both these symmetries are recovered when taking the continuum limit on any reasonable 
discretization scheme. And thus, an accurate enough partition must yield suitable results. 

The stationary probability distribution for the KPZ problem in ID is known to be 

m 

^stat[/i] ~ exp 1^ y dx {d^hf^ . 
For the simplest discretization scheme in Eq. f l22|) . we have 

Vs^^m ~ exp 1^ E ^ [iHrr + i^n'] | • (52) 

Inserting this expression into the stationary Fokker-Planck equation several terms cancel, 
and the ones surviving can be expressed as 

^-j:l\iHt'r+iHry]Lf. (53) 

Clearly, the continuous limit of this expression is of the form 

\v J dx (dxh)"^ d^h, 

15 



that, as is well known is identically zero. A numerical analysis of Eq. (153|) indicates that 
this expression is several orders of magnitude smaller than the value of the pdf 's exponent 
[Eq. and typically behaves as 0{1/N), where N is the number of spatial points used 

in the discretization. Moreover, using expressions with higher accuracy for the differential 
operators one gets an even faster approach to zero. This indicates that the problem with the 



fluctuation-dissipation theorem in 1 + 1, discussed in [ll|, ll5| can be just circumvented using 
more accurate expressions. It is also worth commenting that, if a consistent discrete scheme 
is built from the discrete scheme in [ll|, it would also violate the fluctuation-dissipation 
relation. 



VI. ON THE VARIATIONAL FORMULATION OF KPZ 



In Ref. 26|, a variational formulation was introduced for the KPZ equation. There it was 
shown that Eq. ([T]) can be written as 



where (for F = 0) 



dn[x, t) 
A2 



g[h] = |^e^"(^'*)— 
and the function T{h) is given by 



(54) 



(55) 



m 



X 



The way in which the functionals J^[(f)] and Q[h] are related is also shown in Ref. 26|. It is 
also easy to prove that the functional Q[h] fulfills the Lyapunov property dtQ[h] < 0. 
According to the previous results, we can write the discrete version of Eq. f l55|) as 



^ ^ 8z/Y 2 



Now, introducing this expression into the discrete version of Eq. 0541) . and through a simple 
algebra, we reobtain Eq. ([7]). This reinforces our result, and clearly indicates the need to be 
consistent when considering a discrete version of the KPZ equation. 
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VII. SOME NUMERICAL RESULTS 



We present here some results obtained by numerically integrating the KPZ equation in 
ID. Our aim is to compare the standard discretization scheme [Eq. flTTj) ] with the consistent 
ones presented in Eqs. ( ITU]) and ( 15^ . 

To solve Eq. ([1]) we discretize t) along the substrate direction x with lattice spacing 
a = 1. We employ a second-order Runge-Kutta algorithm (see e.g. [34]) with periodic 
boundary conditions. Then the equation of motion 

h, = V L{h,) + \^{h,) +50 = Pihj) + ^ i3 (56) 

is integrated according to the recursive relation 

hj{t + At) = + ^i9i+ 92) + {^tf'^Uj, (57) 

with 

91 = F{hj) 



g2 = F{h,+Atgi + {AtY/^ + uj), 

where Uj is a Gaussian random variable. 

Without loss of generality, the interface dynamics can be described in terms of the di- 
mensionless parameter A = {2e/i'^Y^'^\. In practice, we set z/ = e = 1 and allow A to 
vary. 

The numerical results show that the interface scaling does not depend on the discretization 
scheme. As shown in Fig. [1], the dynamics fit into the KPZ universality class for all the 
discretization schemes. The global width, that concerns the fluctuations of the growth 
height around its mean value, scales according to the Family- Vicsek Ansatz as 

Wmt) = t^fit/t,)=t^fit/L^), (58) 

where the scaling function f{u) is defined as 

{const if M ^ 1, ^ ^ 

(59) 
u otherwise. 

On the other hand, correlations can be analyzed in the reciprocal space by means of the 
structure factor 

Sk{t) = (hk{t)h_k{t)), (60) 
17 




FIG. 1: Numerical results for the global width and structure factor (inset) averaged over 100 runs 
in a system of size L = 1024 with A = 4. We employ the different discretization schemes indicated 
in Eqs. ()47p . (08]) and (|30l3ip (from bottom to top). Curves are slightly shifted vertically for clarity. 
Lines are plotted as a guide with exponent /3kpz = 0.3 (dashed), j3 = 0.25 (dot-dashed), and —2 
(solid, in the inset). All discretization schemes are consistent with the KPZ scaling. 

where hk(t) is as before (see Sec. lIV|) the Fourier transform of the interface profile. According 
to the previous scaling Ansatz, S(k,t) scales as -^ith the roughness exponent a. 

We observe that all the discretization schemes are consistent with the KPZ scaling, with 
the KPZ exponents a = 1/2 and /3 = 1/3. It can also be observed from Fig. [T] that 
the crossover from the transient linear (Edwards-Wilkinson) behavior to the asymptotic 
nonlinear (KPZ) behavior appears earlier in both alternative discretization schemes than in 
the standard one [Eq. (jS])]. This effect is presumably related to the fact that the nonlinearity 
of the alternative schemes always makes a much stronger contribution to the dynamics than 
the one in the standard scheme, see Fig. [2] below. This way, the threshold contribution from 
the nonlinearity is received sooner, resulting in an anticipated departure from the transient 
linear regime. 
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FIG. 2: Solid line: Time dependence of the nonlinear contribution in two discretization schemes 
[Eq. (jlOp (left frame) and Eq. (j53p (right frame)]. We also depict the Galilean contribution (circles) 
and the excess of fluctuations (diamonds) as defined e.g. in Eq. (jSip . For both cases, the excess of 
fluctuations are comparable with the Galilean contribution. 



In order to analyze the excess of fluctuations that such discretization schemes present 
with respect to the standard one, we extract the Galilean invariant fluctuations from the 
quadratic term of the equation of motion. In Fig. [2] we depict the time dependence of the 
different nonlinear contributions for both alternative discretization schemes. On the left we 
have the comparison between the discretization scheme Eq. (H7!) and the one in Eq. (l48ll . 
while on the right we compare the scheme in Eq. fHTj) to the one in Eqs. fl30ll3ip . 

The first point to note is the rapid saturation of nonlinearities. This behavior is consistent 
with the rapid saturation of local fluctuations, which behave as G(£, t) ~ t^^ and saturate for 



t > 



35| . We also observe that for both cases, the excess of fluctuations are comparable 



to (or even larger than) the Galilean contribution. As we pointed out in previous sections, 
this result, together with the persistence of KPZ scaling shown in Fig. [H would imply that 
Galilean invariance is not such a fundamental symmetry as usually considered. It is worth 
remarking that the present results are not affected by either changing L or using other 
algorithms 211] to perform the numerical simulations. 



19 



VIII. CONCLUSIONS 



The moral from the present analysis is clear: due to the constraint imposed by the Hopf- 
Cole transformation [Eqs. and (jH])] the discrete forms of the Laplacian and the nonlinear 
KPZ term cannot be chosen independently. 

Although the present work is focused on the relation between the diffusion equation 
with multiplicative noise and the KPZ equation, the consequences of this analysis are more 
general. The discrete versions of any set of related differential equations should be obtained 
taking into account the original (or leading) equation and the transformation rules. It 
is worth remarking here that a related analysis was done in {q], but there the emphasis 
was on the study of the strong coupling limit and the mapping onto the directed polymer 
problem, without commenting at all about the consistency among the discrete versions of 
the differential operators. 

The results discussed here are general; they neither depend on space dimensionality 
nor rely on variational representations. Nonetheless, the recently introduced variational ap- 
proach for KPZ |26| offers an adequate framework in order to make a consistent discretization 
of the KPZ equation. 

Regarding the recently introduced pseudo-spectral approach jl^,[l6|, 36|, in addition to its 
known advantages, in principle, it seems to have the virtue of being "transparent" respect 
to the present problem. In this respect, we have shown the relation that exists between 
it and the present analysis. Besides that, when analyzing inhomogeneous situations where 
defects or impurities are present, such methods do not apply and it is again necessary to 
resort to real-space discrete form of the differential operators j^l, [s^]. Another aspect to 
consider is related to the situation in jssi, where a problem of numerical instabilities (a 
computational problem) in discrete growth models has been tackled by introducing higher 
order contributions (changing the physics of the problem!). It is worth indicating that 
such an instability does not seem to arise (or at least it arises latter) in pseudo-spectral 
treatments of the same problem. Hence, due to the relation among both formalisms, it 
seems reasonable to expect that such instabihties could at least be delayed if a consistent 
discretization scheme, together with higher order discrete operators, is used. 

Regarding the two main symmetries associated with the ID KPZ equation, the 
fluctuation-dissipation relation and Galilean invariance, we have found a couple of relevant 
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results. It is clear that both these symmetries are recovered when taking the continuum 
limit of any reasonable discretization scheme. And thus, an accurate enough partition must 
yield suitable results. 

The fluctuation-dissipation relation essentially tells us that the nonlinearity is not opera- 
tive in ID and for long times. Our analysis indicates that the problem with the fluctuation- 
dissipation theorem in 1 + 1 can be circumvented by improving the numerical accuracy. Or 
this is at least what would happen if the interface were smooth enough. We are not com- 
pletely free of surprises coming from the irregular nature of rough interfaces (as we already 
mentioned we expect a Holder exponent strictly smaller than 1/2 for (i-dimensional KPZ 
interfaces). In any case, our simulations have indicated that our strategy of improving the 
numerical accuracy yields operative results. 

Galilean invariance has been always associated with the exactness of the ID KPZ expo- 
nents, and with a relation that connects the critical exponents in higher dimensions. How- 
ever, it is worth remarking that this interpretation has been recently criticized 2^. Our 
analysis indicates that if the numerical solution obtained with a finite differences scheme 
that is not Galilean invariant yields the well known critical exponents, that would strongly 
suggests that Galilean invariance is not a fundamental symmetry as usually considered. It 
is worth commenting that the results presented here for different consistent discretization 
schemes show all the same critical exponents as the standard one, Eq. f l42|) . 

Here we remark that in the present work we have only emphasized the existing constraints 
introduced by the local transformation on the discrete versions of the differential equations. 
No attempt is made here of choosing the most suitable spatial discretization scheme with 
regard to a given KPZ feature, nor to present a deep analysis of results regarding the violation 
of Galilean invariance. The study of such aspects, together with the evaluation of the effects 
of the relations obtained among the discrete operators on different relevant quantities as 
well as other problems will be the subject of further work. 
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